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Abstract — Traditionally, multi-trial error/erasure decoding of 
Reed-Solomon (RS) codes is based on Bounded Minimum Distance 
(BMD) decoders with an erasure option. Such decoders have 
error/erasure tradeoff factor A = 2, which means that an error 
is twice as expensive as an erasure in terms of the code's 
minimum distance. The Guruswami-Sudan (GS) list decoder can 
be considered as state of the art in algebraic decoding of RS 
codes. Besides an erasure option, it allows to adjust A to values 
in the range 1 < A < 2. Based on previous work |1|, we provide 
formulae which allow to optimally (in terms of residual codeword 
error probability) exploit the erasure option of decoders with 
arbitrary A, if the decoder can be used z > I times. We show 
that BMD decoders with 2:bmd decoding trials can result in lower 
residual codeword error probability than GS decoders with zgs 
trials, if zbmd is only slightly larger than zgs- This is of practical 
interest since BMD decoders generally have lower computational 
complexity than GS decoders. 

I. Introduction 

Multi-trial error/erasure (MTEE) decoding or Generalized 
Minimum Distance (GMD) decoding f2\, (3\ is a technique 
which applies multiple decoding trials of an error/erasure 
decoder on each received word, each time with a different 
number of erased most unreliable symbols. The ideas behind 
this approach are to not let unreliable received symbols inter- 
fere the decoding process and to exhaustively try the set of 
most promising erasure patterns. MTEE decoding performs 
surprisingly well, especially when the channel is in good 
shape. This is naturally the case when we consider concate- 
nated codes. Here, the inner code and the channel can be 
considered jointly as a super channel which, due to the inner 
decoder's error-correcting capabilities, has low symbol error 
probability. 

We investigate a particular concatenated code construction 
which is widely used in practice and standards, e.g. the 
Consultative Committee for Space Data System's (CCSDS) 
Telemetry Channel 14|. In this construction, the inner code 
is a convolutional code with a Maximum Likelihood (ML) 
decoder The outer code is a traditional Reed-Solomon (RS) 
code. We stress that the inner code needs to be tailbited to 
insulate channel error events to single symbols of the outer 
received word. 
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Our target is to minimize the residual codeword error 
probability after decoding. We consider threshold erasing, 
which means that each output of the inner ML decoder is 
measured against a set of z > 1 real thresholds ' , . . . ,Tz , 



(2) 



If the reliability of the symbol (which is 

(^) 



an output of the inner ML decoder) falls below threshold T^. 
in decoding trial k, 1 < k < z, then the symbol is erased 
and replaced by the erasure marker X. The threshold erasing 
method dates back to Blokh and Zyablov |5| and is different 
from the symbol erasing method used in Forney's original 
work about GMD decoding. There, the received symbols 
are ordered according to their reliabilities and an increasing 
number of most unreliable received symbols is erased in each 
of the z decoding trials. 

Currently, the most powerful technique for algebraic decod- 
ing of RS codes is the Guruswami-Sudan (GS) list decoder (|6l . 
It can be parametrized to obtain error/erasure tradeoff factors 
A in the range 1 < A < 2. A expresses the relative cost of 
errors compared to erasures in terms of required Hamming 
distance. Generally, increasing the multiplicity parameter v 
brings along higher list decoding radius, increased decoding 
complexity, and smaller A. We will elaborate the latter fact in 
the course of the paper 

The GS decoder has been extended to a soft-input algorithm 
by Kotter and Vardy in their award-winning 2003 paper 
|7|. Their algorithm is based on setting the multiplicity of 
each interpolation point in the GS decoder according to the 
reliability of the corresponding received symbol. Another 
promising approach for soft-input decoding of RS codes has 
recently been published by Nguyen et. al. [8| and is based 
on rate-distortion theory. In our work, we investigate the 
potential of threshold erasing, when the outer code is decoded 
in multiple trials with the GS decoder The results are based 
on our previous papers |[l], |l9], in which we consider outer 
BMD decoding (A ~ 2) of Bose-Chaudhuri-Hocquenghem 
codes and outer decoding of Interleaved Reed-Solomon (IRS) 
codes (A = + 1)/^, ^ e N \ {0}), respectively. 

For the sake of completeness we should also mention 
other publications on related topics, e.g. maximization of the 
decoding radius of concatenated block codes with an outer A- 
decoder using threshold erasing ifTOl . ifTTl and symbol erasing 
lfT2l - lfT4ll . Outer list decoders have already been considered 



by Nielsen ifTSl . but with the aim of maximizing the decoding 
radius of the concatenated code construction. An overview of 
the different erasing techniques with an arbitrary number of 
decoding trials is given in |16|. 

The rest of the paper is organized as follows. In Section |II] 
we describe structure and threshold-based MTEE decoding of 
the aforementioned concatenated code construction. We use 
results from [T|, f9l to derive optimal threshold locations for 
outer decoding with 1 < A < 2 in Section Hill Here and in 
the rest of the paper, optimal means minimizing the residual 
codeword error probability. Section |IV] deals with the GS 
decoder's non-constant A and shows how our result from 
Section [III] can be applied nevertheless. Optimal threshold 
locations are used in Section |V] to plot error probability curves 
of an exemplary concatenated code. It will turn out, that for 
the considered setting the high-complexity GS decoder is in 
many cases not worth the effort and multiple trials of low- 
complexity BMD decoding yield comparable or even lower 
residual codeword error probabilities. We conclude our paper 
in Section [yl] 

II. MTEE Decoding of Concatenated Codes 

A concatenated code C{n,k,d) consists of an inner code 
C'(F2;n',fc' = TO,rf) and an outer code C°{¥2^^;n°,k°,d°). 
The resulting concatenated code C is binary and, w.l.o.g., we 
restrict ourselves to this most practical case. 

The information vector a° e F^™ is encoded into an 
outer codeword c° := (eg, . . . , c°o_i) £ C° C F^°. of the 
outer code. Each 2™-ary symbol c°, j = 0, . . . ,n° — 1, can 
be interpreted as a binary vector e F^ of length m. 
These vectors serve as information for the inner code and 
are encoded into inner codewords e C C Fj . Arranging 
the as columns of a matrix gives the codeword matrix of 
the concatenated code C, which is transmitted over a binary 
symmetric channel (BSC) with crossover probability p. 

The receiver obtains erroneous columns :— c^ +e^, which 
are fed into the ML decoder for C It returns inner codeword 



estimates c^- := dec'(r^ ). The information parts a^ are ex- 
tracted from the and mapped to symbols r° e F2m . The re- 
sulting vector r° := (rg, . . . , r'^a^i) is the input for the MTEE 
decoder of C°. The MTEE decoder performs erasing with the 
threshold set T ■= \t[''\ . . . ,ri^^|, T^^^ < • • • < Ti^^. It 
calculates a reliability value Vj for every received symbol 
r° = a^ according to 

Pr(r;.|c;-) \ 




.Ee.ec.\{a.}Pnri|c' 
and T is applied in the following manner: 



(^) 



X, 



(^) 



k = 1,...,2;. Note that the particular calculation of the 
reliability value stems from Wl\ Corollary to Theorem 1] and 
results in decision regions which minimize both the error- 



and the error-or-erasure probability of the outer decoder at 
the same time. Result of the erasing procedure is the input 
list I := {??,..., f^}, in which r^=(f^_o,...,f^^„„„i). 
Each element of the input list is fed into the outer decoder 
with 1 < A < 2 and multiplicity v. Since we allow the outer 
decoder to be a list decoder, each decoding trial potentially 
returns a result list pk- These lists are merged into the overall 
result list TZ := lJfc=iPfc- We have a decoding success 
whenever c° G TZ. 

III. Optimal Thresholds Locations 

As a starting point for our derivation of the optimal thresh- 
old locations we generaUze |[T] Theorem 1]. 

Several cases are possible when a single received symbol 
r°, which could be either correct (r° ~ c'j) or erroneous 
(7-° = c° + e°), is considered. First, the symbol might be 
correct and never erased by any threshold. We denote the 
probability of this event by 



Pr 



Pr(r° = c° and never erased) 



Second, the symbol might be erroneous and never erased, the 
probability of this event is 

pi := Pr(r° ^ c° and never erased). 

Third, the symbol might be erased by every threshold in T, 
in this case we do not distinguish whether it is correct or not 
and denote the probability by 

Pc := Pr(r° always erased). 

The last two cases are for correct and erroneous symbols 
that are not erased by thresholds T^^\ . . . ,Tj,^\ but erased 

(z) (z) 

by all (larger) thresholds T^^'j^ , . . . , '. The corresponding 
probabilities are 

:= Pr (r° = c° and erased by r^^[_\ but not by xj^"'' 

p^ := Pr (^r° c° and erased by r^^^\ but not by T^^''> 

It is clear that these probabilities must sum up to one, i.e. 

Pr+Pl+Pc + Efc=l(Pfc +£fe) = 1- 

Since it is similar to the derivation of [1, Theorem 1], we 
omit the generalized derivation here and immediately state the 
following theorem. 



Theorem 1 If the outer decoder has error/erasure trade- 
off factor A, 1 < A < 2, and can correct up to ( in- 
cluding) 6 erasures, then the following conditions are nec- 
essary and sufficient for an optimal MTEE threshold set 



Pi 

Pc 



Pc, 



[prvif--^. 



and 



Wk = 1, 



Pk" Pk 



^ Ek+l Pk+l- 



For T fulfilling these conditions, the residual codeword error 
probability can be approximated by 



=p^^iP^Pl)'^'- 



^ipf'p^y^'-^^- (1) 



In case of BMD- and many other decoders S = d° — 1. 
However, we will see later that for the GS decoder we might 
also require smaller values of S. 

Following |5|, we state simple approximations for the 
probabilities of Theorem [T] in our previous paper [fl]. We 
repeat them in Lemma [T] to clarify the further derivation of 
the optimal threshold set. The lemma is based on spherical 
approximations of the inner code's Voronoi cells and the 
exponential error bounds for erasure schemes derived by 
Forney ifTTl . which generalize Gallager's error exponents for 
the BSC nil. 

Lemma 1 (Senger et. al. lU) Simple approximations of the 
probabilities pc,pi,Pf., and p^ are given by 



Pc 

Pi 
Pk 



Ek 



exp(-(i;o(i?') 
exp(-(i;o(i?') 
exp (- (So(i?') 
exp(-(So(i?') 



s Ti^^' 1 n- 



sT, 



k+l 



s tI^'^ 1 n- 



fc = 1, . . . , z — 1, where Eq{R^) is Gallager's error exponent 
for ML decoding of a code with rate i?' and transmission over 
a BSC. s, < s < 1/2, is the corresponding optimization 
parameter 

The conditions from Theorem [T] and the approximations 
from Lemma [T] allow to obtain analytic formulae for the 
optimal threshold locations. Their number z, the rate i?' of the 
inner code and A are parameters. Inserting the approximations 
into the conditions results in the following system of z 
recurrent equations. 



Pt = Pc 



(A + l)n 



(A 



(2) 



(3) 



and, Vfc = 1, 



2, 



Pk'" Pk 



Pk+l Pk+l 



A- 1 



k+2- 



(4) 



Equations (|2]), (O, and ^ allow to prove our main theorem. 



Theorem 2 The optimal threshold set T — \ 



for MTEE decoding of a concatenated code with an inner 



ML decoder and an outer decoder with error/erasure tradeoff 
factor A, 1 < A < 2, is given by 

where Eq{R'-) is Gallager's error exponent for the BSC, s is 
the corresponding optimization parameter, < s < ^, and 



F(A) 



A-l 



- A 



\ z-l 

2(1^ -A 



(5) 



Proof: The statement follows from the unique solution of 
the recurrence relation (|2), (l3), and (|4) for 1 < A < 2. ■ 

Corollary 1 For outer BMD decoding, i.e. A = 2, the optimal 
threshold set is given by 



T, 



(,) ^ go(ff)(2fc-l) 
s(2z + l) 



Proof: The statement follows from the unique solution of 
the recurrence relation ([2), ©, and (|4) for A = 2. ■ 

Corollary [T] coincides with a result of Blokh and Zyablov 
||5]. Thus, we obtain their threshold location formula as a 
special case of our main Theorem |2] 

Fig. [T] shows the optimal threshold sets for z — 20, i?' — 
1/2, p = 0.02, and A = 1.1, . . . , 1.9, 2.0. Each line represents 
one threshold set. Darker color of the curve means larger A. 
The optimal threshold set for outer BMD decoding (A = 2, 
see Corollary [TJ is given as a reference. Note that F{X) is 
constant for fixed A and z, other crossover probabilities p of 
the BSC simply scale the threshold locations by a factor 




20 fc 



Fig. 1. Optimal threshold sets according to Theorein|2]and Corollary [T] for 
z = 20, = 1/2, p = 0.02, and A = 1.1, . . . , 1.9, 2.0. 

It is easy to prove that T^^' is non-increasing with de- 
creasing A, a fact which can also be observed in Fig. [T] This 
means that with decreasing A, the number of erased symbols 
generally becomes smaller We could have expected such a 
behavior since with decreasing A, the relative cost of errors 
decreases and thereby also the effect of erasing unreliable 
received symbols. 

We can use Theorem [T] to obtain an approximation of the 
residual codeword error probability after MTEE decoding with 



an optimal threshold set obtained by Theorem |2] To do so, we 
use the second term from ([T) and write 

Inserting the approximation of pi from Lemma [T] gives 



Pe- (exp(-(£;o(ff) + sT, 



in which we can replace tJ^'' as given by Theorem |2] for 
1 < A < 2 or Corollary [T] for A = 2, respectively. This results 
in the following theorem and its corollary. 

Theorem 3 The residual codeword error probability ofMTEE 
decoding of a concatenated code with an inner ML decoder, an 
outer decoder with error/erasure tradeoff factor A, 1 < A < 2, 
maximal number of correctable erasures 5, and an optimal 
threshold set f = ^[^'\ . . . , ri^-*! can be approximated by 



Pe ^exp -2Eo{R')5 



- 1 



-A 



Corollary 2 For traditional outer BMD decoding, i.e. \ — 2 
and d = d° — 1, we have the approximation 



Fe»exp -2£;o(i?')(d°-l) 



2z + l 



So far, we assumed that A is constant for any number of 
erased symbols. This is true for BMD decoders but not for the 
GS decoder as we will see in the following section 

IV. Dealing with the GS List Decoder's 
Non-Constant A 

The decoder capability function {DCF, a constraint on the 
number r of erasures and the number e of errors, that can be 
corrected concurrently) of a BMD decoder is 

n° -T-2e>k° ^ 1. 

For T erasures, < t < d° — 1, the decoder fails to correct 
£bmd(''') '■— {n° — fc° + 1 — t)/2 or more errors. The indeed 
constant A for any number of erasures is given by the negative 
reciprocal value of eBMD(T)'s slope, i.e. 



A 



BMD • — 



de 



BMD 



dr 



= 2. 



The situation is different for the GS decoder For simplicity, 
we restrict ourselves to the best (in terms of achievable list 
decoding radius) case, i.e. multiplicity u ^ oo. It's DCF is 

{n — T — e)^ 



resulting in EgsIt) — 
Ags(t) := - 



> fc - 1, 
^(fc° - l)(n° 



and 



f decsit) 



V dt 



1 



2V(fc-l)(n-r) 



which is a strictly monotonic increasing function of r and 
thereby not usable in Theorem |2] We will now show that 
near-optimal threshold locations for the GS decoder can be 
calculated using Theorem |2] 

It is straightforward to see that for any r, a decoder with 
radius eGs(''') can be transformed into a decoder with radius 
^Gs(^) £gs{t) by simply discarding all decoding results 
with T erasures and e > £gs(^) errors. This fact and the 
monotonicity of A(r) allow to conclude that any tangent of 
£gs(''") at t = k, < k < (i° — 1, specifies a tangent decoder 
with radius 

deGs(t) 



dt 



(r-«) 



and constant error/erasure tradeoff factor 

AgS,k 



dt 



t — K. 



= 1- 



k° - 1 



2^(fc° - l)(n° - k) J 

that can be imitated by the GS list decoder. Its maximum 
number of correctable erasures (5gs,k is obtained by solving 
^gs,k{t) — for r and taking the floor of the result. 

Since Ags.k is independent of t. Theorems |2] and [3] can 
be applied with Ags,k and (5gs,k to obtain optimal threshold 
locations and residual codeword error probabilities for tangent 
decoders which can be imitated by the GS decoder The 
optimal tangent decoder is determined by 



arg mm 

0<K<d°-l 



3GS,k ■ 



VAgs,„-1 J 



(ags^-i) 



-A 



(6) 

which is independent of the the ML error exponent. Thus, 
tangent decoders determined by (|6]l are optimal for all BSC 
crossover probabilities. 

V. Simulation Results - Traditional BMD 
Decoding can Beat the GS List Decoder 

Let us consider the outer RS code C°(F28; 255, 144, 112) 
with GS decoder We consider zgs '■— 1, 5, 10 outer decoding 
trials. Based on (|6]l. Table J] states the parameters of the 
corresponding optimal tangent decoders. 



ZGS 


K* 




<5gs,k* 


1 


41 


1.69126 


107 


5 


72 


1.79208 


110 


10 


85 


1.84699 


111 



TABLE I 

Optimal tangent decoders for zqs 



1,5, 10. 



As inner code, we assume a tailbited rate i?' = 1/2 
convolutional code with ML decoder. This allows to use 
Theorem |3] in order to plot the solid red residual codeword 
error curves for outer GS decoding in Fig. Additionally, we 
consider outer BMD decoding and allow the decoder to be run 
^BMD ■= 1, 5, 10 times (dashed blue curves). We observe that 




Fig. 2. Residual codeword error probability curves for C°(F2s;255, 144, 112) and i?' = 1/2. 



the gain of tangent decoding diminishes for growing z. Since 
both residual codeword error probabiHties (optimal tangent 
decoder and BMD decoder) converge to the same value, i.e. 

Pe'^exp(-£;o(ff)(d°-l)n') 

(dash-dotted black curve), we conclude that for every number 
ZGS of outer GS decoding trials, there exists a number zbmd > 
2;gs of outer BMD decoding trials that achieves either the same 
or lower residual codeword error probability. This allows to 
trade a number of high-complexity GS decoding trials for a 
(generally larger) number of low-complexity BMD decoding 
trials, extending the options of the system designer. 

VI. Conclusions 

We generalized our results from |IT] to the case of arbitrary 
error/erasure tradeoff factors A in the range 1 < A < 2. 
We derived formulae for optimal thresholds applicable in 
MTEE decoding, our generalization allows to use the GS 
list decoder for the outer code. Based on our derivation, we 
gave approximations of the residual codeword error probability 
after outer decoding for the full range of A. This allowed 
to compare outer GS list decoding with traditional, low- 
complexity, BMD decoding. Our main result is that for the 
particular concatenated coding scheme under consideration 
(outer RS code, inner convolutional code with ML decoding, 
e.g. used in pl ), zbmd trials of outer BMD decoding can 
outperform zqs trials of GS decoding if zbmd > ^gs- This 
is interesting for practical applications, since BMD decoders 
have low computational complexity and are widely deployed. 
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